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It is shown that the quasi-universal ratio q — liruy^o & S/ C ~ ±1 of the Seebeck coefficient 
to the specific heat in the limit of T = observed in a series of strongly correlated metals 
can be understood on the basis of the Fermi liquid theory description. In deriving this 
result, it is crucial that a relevant scattering arises from impurities, but not from the mutual 
scattering of quasiparticles. The systematics of the sign of q is shown to reflect the sign 
of the logarithmic derivative of the density of states and the inverse mass tensor of the 
quasiparticles, explaining the systematics of experiments. In particular, the positive sign of 
q for Ce-based and / 3 -based heavy fermions, and the negative sign for Yb-based and f 2 - 
based heavy fermions, are explained. The case of non-Fermi liquid near the quantum critical 
point (QCP) is briefly mentioned, showing that the ratio q decreases considerably toward 
antifcrromagnetic QCP while it remains essentially unchanged for the ferromagnetic QCP 
or QCP due to a local criticality. 

KEYWORDS: Seebeck coefficient, specific heat, Fermi liquid, strongly correlated metals, peri- 
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Strongly correlated electron systems exhibit some quasi-universal ratios among various 
physical quantities. The Wilson ratio, 1 xlli °f the magnetic susceptibility \ to the Sommer- 
feld coefficient 7, and the so-called Kadowaki- Woods ratio 2 A/7 2 , A being the coefficient of 
the T 2 -term of the resistivity, take a quasi-universal value in a series of strongly correlated 
electron systems up to logarithmic accuracy. The validity of the Fermi liquid theory and 
the dynamical nature of the mass enhancement lie behind the universality. 3 ' 4 Quite recently, 
Behnia, Jaccard, and Flouquet have revealed that the ratio of the Seebeck coefficient S to 
the specific heat 7T takes a quasi-universal value, in the low-temperature limit, for a series 
of correlated compounds. 5 The purpose of this Letter is to clarify the reason for the universal 
value for eS/'yT on the basis of the Fermi liquid theory for thermoelectricity. 6,7 

A firm starting point is that the canonical formula for the Seebeck coefficient S in the zero- 
temperature limit 8 is valid even though strong correlations are apparent among the electrons: 



1/12 



J. Phys. Soc. Jpn. Letter 

where e (> 0) is the elementary charge, a(e) is the conductivity of electrons with energy e, 
and /j, is the chemical potential at T = 0. Although the above expression was derived for a 
single-band Fermi liquid, 6 ' 7 there is no reason for doubting its validity also in the multiband 
Fermi liquid. 9 Therefore, we discuss the Sccbcck coefficient on the basis of formula (1). Then, 
we use the conductivity for a multicarrier system: 

a(e) = J>(e), (2) 

i 

in the expression of the Seebeck coefficient, (1). Here, <jj denotes the conductivity due to 
the i-th Fermi surface of multicomponent carriers. By some algebra, we obtain the rule of 
summing the Seebeck coefficient of each component as follows: 

a 

Single-Band Fermi liquid 



First, we discuss the case of heavy fermions with a single Fermi surface. The conductivity 
of a single-band Fermi liquid is approximated as 4 

a(e) = e 2 r(e) £ 5(e - e(k))(v a (k)v* a (k)), (4) 

k 

where v* denotes the renormalized velocity by vertex corrections due to the many-body effect, 
and (• • • ) denotes the average over the direction a (= x,y,z). Here, v(k) = dE^/dk, E? 
being the dispersion of the quasiparticle, and r(e) denotes the transport collision time of the 
quasiparticles with energy e. Hereafter, the renormalized velocity v& is approximated by v a , 
because the impurity scattering plays a dominant role for which such renormalization is not 
crucial. 

The transport collision time for the Fermi liquid is given, around the Fermi level in the 
zero-temperature limit, as 

r(e)- 1 =r imp (e)- 1 + B(e- f i) 2 , (5) 

where Ti mp (e) is the collision time due to impurity scattering, and B is a constant of the 
order of the inverse effective Fermi energy. The second term in the collision rate (5) makes no 
contribution to S since the following relation holds: 

lim — lnr(e) = lim — lnr imp (e). (6) 
oe e-^tt oe 

Therefore, the Seebeck coefficient for a single-band Fermi liquid is expressed as 

I 



Tr 2 klT,. Id, 
S = - — -^-lim — lnr imp (e 
3 e e=M oe 



Zs&(t-E R ){dv a (k)/dk a ) 
ZkS(e-Er)(v a (k)v a (k)) 
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^—T. 



(7) 



The collision rate due to impurity scattering in the i-matrix approximation is given by 

2irz 2 {e)N*{e)u 2 



imp 



x c 



l + [7rz(e)N*(e)u] 2 '^ imp ' (8) 
where Cj mp is the impurity concentration, N*(e) is the density of states (DOS) of the renor- 
malized quasiparticles, and z(e) is the renormalization amplitude. 10 In deriving eq.(8), we 
have assumed the s-wave impurity potential and also neglected the real part of the electron 
self-energy due to impurity scattering. Then, the logarithmic derivative of T imp is given as 

M/x)7v*( M n 2 



lim — In t; 



imp 



(0 



1 x |-ln7V*(e) 



(9) 



[7rz(fi)N*(fx)u} 2 + l de 
Here, we have used the estimation lim e ^ \dlnz(e)/de\ <C lim^^ \dhi N*(e)/de\, which has 
been verified by explicit calculations, e.g., FLEX approximation for t-t' Hubbard model near 
the half-filling. 11,12 It is noted that z(/i)iV*(//) is of the same order as the unrenormalized one 
since the mass enhancement in N* is cancelled by the smallness of z, while the logarithmic 
derivative of N* is highly enhanced in a heavy fermion situation. 

In order to make the discussion explicit, we proceed with the periodic Anderson model 
(PAM): 



I' 



(10) 



where the notations are conventional. According to the Fermi liquid theory based on PAM, 4 
the effective Hamiltonian of quasiparticles near the Fermi level is given as 



H, 



qp 



E 



erc+ c z + i f ft fr + \Vtct ft + h.c.l } , 



(11) 



where the renormalized /-level is determined by the /-electron self-energy Ef and the renor- 
malization amplitude af as, 



i f = a f [e f + E f (n)}, 



(12) 



and Vt = ^/afVr. The fc-dependence of Ef should be rather weaker than \V\ , because 

it is the only way to realize the heavy fermion state by satisfying the two conditions, the 
Landau-Luttinger sum rule 13 ' 14 and the requirement that the /-electron number rif at each 
site be nearly unity, rif ~ 1. Then, the dispersion of the quasiparticles is given as 



E k= 2 



e? + £f 



(13) 



where we have assumed the Fermi level is located at the bonding band which is considered to 
correspond to the Ce-based heavy fermions. 
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Then, the DOS of the quasiparticles, N*(e) = 2^g<5(e — E^), is given by the standard 
calculation as follows: 

N*(e) = 2N cmd € *^~ 2e + ii , (14) 
£{- e 

where AT cond i s the DOS of the conduction electrons, and e*(e) is determined by the following 
equation: 

1 



e*(e) + h ~ J(e*(e)-it) 2 + MV n \ 



(15) 



and is of the order of D, half the bandwidth of conduction electrons. Since | e* (e) — n\ » \i{ — fj,\, 
the DOS is given near the Fermi level, e = fi, as 

N*(e) ~ 2N cond ._ D , = - 6 , , (16) 
\i{ — e\ \i{ — e\ 

where b is a constant of 0(1) because N cond D ~ 1. Thus, the logarithmic derivative of DOS 
is given by the following simple expression: 

1 



J^lniV*(e) 



e=M M - £ f 



(17) 



It is noted that the denominator is of the order of the effective Fermi energy. 

Now we discuss the second term in the brace of (7). With the use of the dispersion relation 
(13), the velocity of quasiparticles v a (k) = dE^/dk a is estimated near the Fermi level as 

a[ ' it + H - 2E % dk a ~ e*00 8k a ' U8j 
where we have considered the fact |ef — E^\ <C e*(/x). The inverse mass tensor is given as 

dv a {k) ^ i t -Eg (desV it-Ez &ej 



(20) 



3/c a e*( / u) 2 \dk a J e*(fi) dk 2 . 

Therefore, the second term in the brace of (7) is calculated as 

J2r6(e-E R )dv a (k)/dk a 2-r, 



ZkS(e-E R )v a (k)v a (k) e-Sf 

where rj = e* (fi)(d 2 e^/ dk^) / (de^/ dk a ) 2 , which is equal to 1/2 for the free electron dispersion 

and for the linear dispersion of conduction electrons. 

Collecting the relations (9), (17), and (20), the factor T in eq.(7) is expressed as 

?,- V + {l- V )[iTz(ix)N*(pi)u] 2 1 

1 + [irz{fi)N*{fi)u] 2 n-if' 1 ' 

Therefore, the ratio q of the Seebeck coefficient given by eq.(7) to the specific heat 7T is 
reduced to the concise form of 

= S e = 3-r ? + (l- ?? )M / u)iV*( / uH 2 I 
q ~Tj 1 + [Kz(fi)N*(n)u] 2 6' 1 ' 

where 7 = ir 2 k 2 i N*(iJ,)/3 is the Sommerfeld coefficient, and b = N*(fi)(if — fj,). In deriving 
eq.(22), we have used eq.(16). The ratio q depends on the character of the impurity scattering: 
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For a weak impurity potential producing the Born scattering {zN*u <C 1), q — (3 — rj)/b, 
while q ~ (1 — n)/b for the impurity with the unitarity scattering (zN*u S> 1). The strength 
of scattering by actual impurities extends from weak to strong coupling. Therefore, we should 
have taken an average over the impurities when we derived relation (8), rather than simply 
multiplying the impurity concentration Ci mp . If we do so, however, a simple expression as (22) 
will not be obtained. Therefore, for simplicity, we simply take the average of expression (22) 
over the strength of the impurity potential. As expected physically, most impurities will have 
an intermediate character. Therefore, it is reasonable to expect that the ratio is of the order of 
0(1), leading to the quasi- universal value for q ~ +1 for Ce-based heavy fermions. A crucial 
point here is that strong renormalization effects in S and 7 cancel each other out, as in the 
case of the Kadowaki- Woods ratio. 3,4 

In the case of Yb-based heavy fermions, we should apply the hole version of the above 
discussions, in which the Fermi level is located on the antibonding band in an electron picture. 
Even in this case, expressions (14), and (16)-(21), are valid. Since > if, the sign of (22) 
changes, leading to q ~ —1. These results explain the observed universal behavior of the 
Yb-based heavy fermion compounds. 5 

Non-Fermi Liquid near Quantum Critical Point 

In the case near the antiferromagnetic (AF) quantum critical point (QCP), where the 
quasiparticles are still well defined, the DOS N*(e) is nearly symmetric around the chemical 
potential. Indeed, near the AF-QCP in 3d, N*(e) oc {const. - [(e - fi) 2 + u$ 1/4 } x N{ oc (e), 
with ujq being the energy scale inversely proportional to the staggered susceptibility and 
-^ioc( e ) being the DOS renormalized by a local correlation effect as discussed above. It is 
noted that ujq has T-dependence, oc T 3 / 2 , in general. Therefore, lim^^ d In N* (e) /de = 
lim e ^ M 31n N\ oc (e)/de. Furthermore, the second term in eq.(7) is shown to remain unenhanced 
by the AF critical fluctuations. Indeed, the velocity v a appearing in eq.(20) depends crucially 
on the position on the Fermi surface, and vanishes in proportion to the inverse of the loga- 
rithm of its energy on the hot lines but remains finite otherwise. Then, the summation with 
respect to k on the Fermi surface remains of the same order of magnitude as away from the 
QCP. 15 

Therefore, the ratio q near the AF-QCP is given as \q\ ~ N* oc (e)/N*(e)\ e=ll , which is far 
less than unity at AF-QCP because N*(/j.) is enhanced compared to N\ oc (fi) by AF critical 
fluctuations, although it does not diverge. On the other hand, near the ferromagnetic (F) 
QCP, q remains of the same order as that away from QCP, because the second term of eq.(7) 
is proportional to the mass enhancement factor both due to local correlation, given by eq.(20), 
and the F critical fluctuations which enhance the effective mass equally at all points on the 
Fermi surface, in contrast to the AF critical fluctuations discussed above. The estimation of 
the first term in eq.(7) is also valid for F-QCP. Then, the ratio q exhibits only a small decrease 



5/12 



J. Phys. Soc. Jpn. Letter 

leaving q ~ ±1 with a logarithmic accuracy. 

Thus, the ratio q considerably decreases toward the AF-QCP but does not change appre- 
ciably for the F-QCP. The estimation for F-QCP is also valid for the local quantum criticality, 
where all the points on the Fermi surface are subject to the effect of critical fluctuations. 16 ' 17 
These predictions may be explored by experiments around QCP tuned by altering the pressure 
or the magnetic field. 



Multiband Fermi liquid 

Actual heavy fermion compounds have a multiband structure and plural Fermi surfaces. 
Therefore, one should calculate Si, the contribution from the i-th band, and sum up by using 
formula (3) with weight Oifo. In the zero-temperature limit, where only the impurity scattering 
is relevant to the thermoelectricity, ai/a is of the same order for all of the bands because the 
renormalization effect cancels out for the impurity scattering 18 unless a renormalization of 
the impurity potential due to the quantum critical fluctuations develops. 19,20 Indeed, the 
conductivity of light (£) and heavy (h) bands is given by 

* e > h ~ e 2 r imp (4' h ) 2 <' h , (23) 

where the collision time r; mp is given by 

= ^ 2 c^u\z^fN^. (24) 

It is easy to see that the renormalization factors z are cancelled among those included in 
v oc z and Np oc z~ l . Then, a e ~ a h if the impurity potentials for light and heavy carriers are 
comparable as expected. 21 Therefore, the contribution of Si from the heavier band dominates 
because |S h | » \S e \. 

Some heavy fermion compounds such as CeCu2Si2 are compensated metals for which we 
have to take into account the multiband conduction electrons. Nevertheless, enhancement 
occurs in the first and second terms in eq.(7) in proportion to — din N*(e)/de\ e=ll in each 
heavy fermion band, maintaining the validity of the result for a single-band Fermi liquid. 
This is because the Fermi level is located just below the level corresponding to the divergence 
of DOS due to the hybridization effect, which causes the hybridization gap in the case of a 
single conduction band. A crucial point is that the sign of din N*(e)/de and the inverse mass 
tensor do not depend on the sign of the velocity but on the curvature of the dispersion of 
quasiparticles. This explains why the q value is positive for Sr2Ru04 even though the 7-band 
(heaviest band) is electron-like with the negative Hall coefficient, Rn < 0. Indeed, the energy 
derivative of the DOS of the 7-band is positive due to the van Hove singularity located just 
above the Fermi level. 22 It is also the case in (BEDT-TTF) salt reported in ref. 5. Indeed, 
the DOS in the Hiickel type tight-binding theory seems to explain the difference in the sign 
of the Seebeck coefficient between the b- and c-directions. 23 
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A prediction based on the above result is that CeRu2Si2 with q ~ +1 without the magnetic 
field H = changes the sign of q just above the metamagnetic field Hm for the temperature 
gradient along the c-axis. This is because the singular peak of the DOS arising from the flat 
band structure around (0, 0,7r/c), which is located above the Fermi level, 24 ' 25 is expected to 
shift down to the Fermi level at H > Hm, leading to the reversal of the signs of din N*(e)/de 
and the inverse mass tensor. 



f 2 -Based Heavy Fermions 

The above theory for heavy fermions is valid for / 1 -based compounds. On the other hand, 
/ 2 -based heavy fermions, such as U- and Pr-based compounds, have a different structure of 
the quasiparticle band. In such systems, one has to take into account the plural /-orbitals split 
by the crystalline electric field (CEF) effect. It gradually became apparent in the mid 90s that 
the effective CEF splitting is considerably suppressed to be less than the renormalized energy 
scale, the effective Fermi energy, while that of the / 1 -based systems is slightly enhanced due 
to the correlation effect. 26,2 Indeed, this is the only way for the mass of quasiparticles to be 
highly enhanced by satisfying the two conditions, the Landau-Luttinger sum rule 13 ' 14 and the 
requirement that the /-electron number na of the i-th orbital (i = 1, 2) at each site be nearly 
unity, rijf ~ 1, as seen below. It was shown on the basis of the investigation of the f°-f l -f 2 
model 27 that the mass enhancement arises only if the /-electron number is nearly 1 or 2, as 
long as the two /-orbitals with a low- lying CEF level are relevant. 

The effective Hamiltonian for the quasiparticles of / 2 -based heavy fermions is given in the 
form: 

k y * =1 > 2 

+ E feu^+Hl (25) 

i=l,2 ) 

where the renormalized /-levels ia's and the renormalized hybridizations V^'s are given by 
the orbital-dependent /-electron self-energy and the renormalization amplitude as 

in = a,if[£if + £jf(/i)], (26) 

and = ^/aifV^. The dispersion of the quasiparticles E k is given as three solutions of 

J] (i* - E k ){e % -E k )=Y J ~ E k ). (27) 

i=l,2 i=l,2 

The schematic behavior of the dispersion is shown in Fig. 1. The renormalized /-levels are lifted 
from their original positions in such a way that the total electrons per site are distributed to 
two /-electrons and the rest to the conduction electrons. The number of correlated electrons 
is counted using the Landau-Luttinger sum rule, while that of the conduction electrons is 
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counted by summing up the fc-points below the Fermi energy. Then, it is easily seen that the 
two /-electrons occupy the orbital of the lower level in unless the renormalized level splitting 
A = £2f — in is far less than the width of each band ~ af\V\ 2 /D. On the other hand, if A is 
significantly smaller than a,f\V\ 2 /D, it is possible to satisfy the above two conditions because 
the lower level then consists of both /-orbitals with nearly equal weight. This is the reason why 
the renormalized level splitting is considerably reduced to make the / 2 -based heavy fermion 
band. 26 ' 2 It is not difficult to see that eq.(27) has a trivial solution at some k as 

E k = e k = e ^^ f (28) 

where V£ = (V? ±V?). 




Fig. 1. Schematic view of the dispersion of / 2 -based heavy fermions in some direction in the Brillouin 
zone. The dispersion of conduction electrons is assumed to be linear, and k-p and fee denote the 
Fermi wavevector and the zone boundary, respectively. The dashed curves are for a noninteracting 
system, and the solid curves are for a renormalized dispersion. The energy scale of the right-hand 
panel is enlarged compared to the left-hand one. It is noted that effective CEF splitting is reduced 
to be smaller than the characteristic energy scale Tq. 



8/12 



J. Phys. Soc. Jpn. Letter 
Formal solutions of eq.(27) near the renormalized /-levels are given similarly to eq.(13) as 

V 2 



E 



± 



e + 



2V 2 



-A- 



V 



2e k \ 4 



A 2 



+ 



2e k 



- —A 

2e k 



(29) 



where e k = e k — E k and \e k \ ~ D except for the narrow region \e k — e\ < max(A, V 2 /D). In 
the situation shown in Fig. 1, the Fermi level is located in the band E~ with e kF < e, which 
simplifies the following analysis. In order to avoid the double occupancy of the orbital- 1 with 
a lower CEF energy, e kp , the energy of the conduction band corresponding to the Fermi wave 
number should be slightly lower than e so that |e& p | » |ef — E^\, |V+|, because the weight in 
the lower level en of the orbital-1 is slightly larger than that of orbital-2. 
The DOS near the Fermi level is given as 

Vl 



N*(e) ~ 4iV 



cond 



' vl 



£+ 2 + 2(e k -E~) 



e — e 



(30) 



where A = AVl/V 2 , and the term of 0(V 2 /(e k - E k ) 2 ) has been neglected. Then, the 
logarithmic derivative of DOS at the Fermi level is expressed as 



-lniV* (e) 



(e - /x) 2 + (e - n)(T - A ) + ](T - A ) 



To 



-A, 



T 



(i-fj,)(i-H + 



T -A c 



4J 



(31) 



where the characteristric energy scale To is defined as To = V+/|efc F |. The inverse mass tensor 
at the Fermi level is given as 



d 2 E k 



dk 2 



x 



I 1 



1 + 



2 ' N . T - A 
£ - A* H o 



T (T - A 



A r 



(32) 



With the use of expression (32) and a similar expression for the velocity dE k /dk a , the second 
term in the brace of eq.(7) is expressd as 

(~ , T -A V 

Y.%^-Efi)dv a {k)/dk 



Y, k 6(e - E^)v a (k)v a (k) 



e-^ + ^^Y + luTo-Ao) 



e=/j. 



To 



(e- fi)[e- fi + 



A r 



AH 
Th 



> 0. 



(33) 

Thus, with the use of expressions (31) and (33) , the factor T in eq.(7) is expressed as 

> — — — — '^-7—^— — — ■ >0 (34) 



4 H + J+{H- J)[irz(j,)N*^)u] 2 
To l + [7Tz(fi)N*(fi)u} 2 * • 



Since H and J are positive constants of 0(1) and N*{^x)To ~ 1, the quasi-universal value of 
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q ~ — 1 follows as in the argument for the single-band Fermi liquid. This result explains the 
fact that q ~ — 1 for UBei3 and URu2Si2- 5 A straightforward prediction is that the Pr-based 
heavy fermions with the filled Skutterudite structure, such as PrFe4Pi2 and PrOs4Sbi2, will 
exhibit the universal ratio q ~ — l. 28 

In contrast, UPd2Ai3 with q ~ +1 is considered to be in the /^configuration 29 and the 
situation would be the same as in the Ce-based compounds, which is consistent with q ~ +1. 
This is because three /-electrons exhibit itinerant-localized dual behavior: two of them are 
localized to form an / 2 CEF level structure and one of them forms the quasiparticles. 30 
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lating discussions. One of us (K.M.) benefitted from a lively conversation with D. Jaccard and 
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